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Abstract 

We investigate the impact of recent direct supersymmetry (SUSY) searches on a non-universal 
sfermion mass scenario focusing on naturalness. One of the advantages of this scenario is that the 
non-universality between third generation and first two generation sfermion masses can relax the 
tension between naturalness and constraints from flavour and CP violating observables. In the 
parameter region, where various phenomenological constraints are satisfied, the constraints to 
this scenario from ATLAS 165 pb _1 "O-lepton" search and 35 pb _1 "6- jet" search are much weaker 
than those to the constrained minimal SUSY standard model, because of differences in the main 
SUSY production processes and the main decay chains. Naturalness can be easily achieved in 
this scenario in accord with the current direct SUSY searches. An additional dedicated analysis 
may be needed to discover /exclude this scenario. 



1 Introduction 



Weak scale supersymmetry (SUSY) is one of the most promising candidates for physics beyond the 
Standard Model (SM). The minimal SUSY extension of the Standard Model (MSSM) has several 
attractive features. For instance, it improves the gauge coupling unification indicating grand unified 
theories (GUT) and provides a dark matter candidate as the lightest supersymmetric particle (LSP) 
[HO [3]. One of the most attractive and important features of the MSSM is to provide a solution to 
the gauge hierarchy problem. The MSSM removes a quadratically divergent radiative correction to 
the Higgs mass up to the SUSY breaking scale, which is a mass scale of supersymmetric particles. 
An unnatural tuning is therefore not required if this scale is around the weak scale, (9(100) GeV. 
In other words, "naturalness" of the theory disfavours the superparticle mass scale ^> (9(100) GeV. 
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However, there may be a tension between experimental constraints and naturalness in the MSSM. 
Flavour Changing Neutral Current (FCNC) processes and CP violating (CPV) observables place 
stringent constraints on the MSSM parameters. Under these constraints, the mass scale of super- 
particles has to be above (9(100) TeV, otherwise the SUSY breaking parameters must be precisely 
aligned in a flavour and CP conserving manner. 

In addition, new constraints have been reported by the ATLAS and CMS experiments at the 
CERN Large Hadron Collider (LHC). They have conducted direct searches for super symmetric 
particles with = 7 TeV pp collision data. With the 35 pb -1 data recorded in 2010, the strongest 
limit on the constrained MSSM (CMSSM) comes from the ATLAS "0-lepton" search [4J. It excludes 
equal mass squarks and gluinos with masses below 775 GeV at 95% confidence level (CL) in the 
Aq = 0, tan f3 = 3, \i > slice of CMSSM. Recently, a new result has been reported by the ATLAS 
with the 165 pb _1 data recorded in 2011 [5]. It excludes equal mass squarks and gluinos with masses 
below 950 GeV at 95% CL in the CMSSM with A = 0, tan/3 = 10, fi > 0. It may be getting more 
difficult to realise naturalness within the framework of the CMSSM in accord with the experimental 
constraint^]. 

The situation may change if one relaxes the assumption of universality among the scalar fermion 
(sfermion) masses. Strictly speaking, the experiments and naturalness constrain different parame- 
ters. Naturalness indicates the parameters that are sensitive to the low energy up-type Higgs mass 
squared, m^, should not be much larger than the weak scal^E Such parameters are typically masses 
of gluino and third generation squarks. 

In contrast, FCNC and CPV observables provide stringent constraints mainly on the first two 
generation sfermion sectors. The constraints on the third generation sector are, on the other hand, 
not so severe. The introduction of non-universality between the third generation and the first 
two generation sfermion masses is therefore possible without conflicting with the FCNC and CPV 
constraints El EH HH Q2] . 

In the light of this observation, one can assume the first two generation sfermions are much 
heavier than the weak scale. This assumption relaxes the necessity of the alignment mentioned 
above among the SUSY breaking parameters and allows the parameters to have some amount of 
flavour and CP violations. At the same time, the gluino mass and the third generation sfermion 
masses can be assumed to be in the region of the weak scale to keep naturalness. In what follows, we 
call this scenario modified universal sfermion mass (MUSM) scenario. The mass spectrum of MUSM 
scenario is supported by several earlier models O [13l H31 H51 HEl HZ]. Various phenomenological 
constraints, such as FCNC, EDM and colour and charge breaking (CCB), are investigated in Refs 

PSl HS1 1201 12H 1221 - 

It is not clear whether naturalness of this scenario can still be realised under the recent direct 
SUSY search constraints. Those searches have set strong limits on the CMSSM parameter space. 
However, there are remarkable differences between the CMSSM and the MUSM scenario. First, 
the main SUSY production processes are qq, qg and gg in the CMSSM, whilst they are gg and 
t\i\ (&i&*) in MUSM scenario. The main decay chain is also different. It is q — > gj -» Xijjj 
(d ~~ ^ Q3 ~~ ^ Xijj) ^ m q > m g ( m g > m q) m the CMSSM. On the other hand, long cascade chains 
such as g —> t\i —> x[bi -> XiW^bi are anticipated in the MUSM scenario. 

Using recent direct SUSY searches, the aim of this paper is to investigate constraints on the 

1 There are some regions in the CMSSM where naturalness can still be realised without conflicting with the LHC 
constraints [6j [7] . 

2 1 2 2 1 22 

2 A parameter tuning appears in the condition — \™i- m 2 \ = _ m i m 2 _ _ m | _|_ _ m||(9( , \ R ), where 

z 2yl — sin 2 2/3 z tan P 

mi is the low energy down-type Higgs mass squared. The mi is not so sensitive to naturalness due to the 1/ tan 2 /3 
suppression. 
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MUSM parameters that are relevant to naturalness. More precisely, we will set constraints on the 
(7713 — rai/2) parameter plane, where m^ 2 is universal gaugino mass and 7713 is the mass of the third 
generation sfermions involved in the SU(5) 10-plet. We shall use the ATLAS "6-jet" search [23] 
with 35pb _1 data and the ATLAS "0-lepton" search [5j with 165pb _1 data. The former search 
targets large tan/3 scenarios in the CMSSM or cases where gluinos predominantly decay to third 
generation squarks. We shall also use the latter search because it currently provides the strongest 
limit on the CMSSM parameter space. 

The paper proceeds as follows: in Section 2, we clarify the parameters in our analysis and identify 
the region of parameter space where various phenomenological constraints are satisfied. Production 
cross sections of SUSY particles and gluino branching ratios are also discussed. In Section 3, we 
review the properties of ATLAS 6-jet search and 0-lepton search. In Section 4, we describe our setup 
for simulation and formulae for the exclusion p- value. We present the results of the ATLAS 6-jet 
search and 0-lepton search constraints on the MUSM scenario in Section 5. Section 6 summarises 
our results. 

2 Viable parameter space, cross sections and branching ratios 

First, we clarify the parameters in our analysis. As discussed in the Introduction, the MSSM 
parameters can be classified into two groups. One is the group of parameters which are relevant 
to naturalness, and the other is for parameters which are less sensitive to it. The members of the 
first group are ttih u , |/i|, m i R i 171 q 3 an( ^ m i/2> where Q3 is the third generation squark doublet. 
Those parameters are defined at the GUT scale. In this paper, for simplicity, we assume an SU(5) 
GUT [211 ES] relation for the soft SUSY breaking parameters and employ universality for the two 
Higgs soft masses. We then define the following parameters: 

m 3 = mi R rriQ 3 m fR , m H m Hu m Hd . (1) 

The second group involves the trilinear couplings and the other soft scalar masses. We assume 
their universality for simplicity. The same symbols Aq and mo are used for the universal trilinear 
coupling and the universal scalar mass, respectively, as in the CMSSM. 

In summary, our interest region is where the three parameters in the first group are of the order 
of the weak scale and the two parameters in the second group are of the order larger than the weak 
scale. 

Group 1: m^ 2 , 7713, ma — 0(100) GeV 
Group 2: m , (\A \) > 0(100) GeV 

The value of |/i| is determined by the condition of electroweak symmetry breaking. The sign of \i 
remains physical. We fix sign(/i) = +, so that the b 57 constraint is relaxed due to a cancelation 
between the charged Higgs-top quark contribution and the chargino-stop contribution. We do not 
address the b —> 57 constraint explicitly because it depends on the off-diagonal entries of V UR1 
V UL and Vd L in this scenario [TT], where Vj I (/ = u,d and I — L^R) is the unitary matrix that 
diagonalises the Yukawa matrix as Y^ iag = Vj L YfVf R . Finally, we take tan/3 = 10 in our analysis. 
This type of MUSM spectrum can be obtained from an Eq GUT model with SU(2) horizontal 
symmetry [HI El [13JH1- Moderate value for tan/3 (5 < tan/3 < 15) is also a prediction of this model. 
In what follows, we fix 7713 = trh for simplicity and take mo = 1.5 TeV. The value of mo is not 

3 The assumption rriH u — mH d at the GUT scale can be removed in this model. 
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m = 1.5TeV, m Hu = m Hd = m^, m 1/2 = 150GeV, tan/3 = 10, ft > mo = 1.5TeV : = mn d = mz, Aq = -7713 - 350, tan/3 = 10, //, > 

. J . 1 . 1 . 1 . 1 u 400 -1 1 ' 1 ' 1 ' 1 ■ 1 1 1 ' ■ ' — 1- 




m 3 [GeV] m 3 [GeV] 



Figure 1: Left: the excluded regions by the Higgs mass bound (the green region) and by the tachyonic 
stop (the pink region) on the (777,3 — Ao) plane at rai/2 = 150 GeV. Right: viable parameter space on the 
(7773 - m 1/2 ) plane. 



sensitive to our collider study as far as mo > 1.5 TeV because the superparticles with masses 
> 1.5 TeV are almost inaccessible in 165 pb -1 data at — 7 TeV due to the small cross sections. 

Let us identify the region of parameter space where various phenomenological constraints are 
satisfied. There are two major constraints to this scenario. One is a CCB constraint where one of 
the mass squared eigenvalues of scalar tops (stops) ra| is negative. The 7T7q negatively contributes to 
low energy values of m~^ and m~^ via renormalisation group evolution in a two-loop level. Moreover, 
if I At I is large at low energy, one of the eigenvalues of the mass squared matrix for stops becomes 
small. These two effects may drive m~ negative. With fixed mo, the CCB (or tachyonic stop) 
constraint requires large 7713, m^ 2 and small \Aq\. 

The other major constraint is the lower mass bound on the lightest CP-even Higgs, h. In our 
parameter region, the lightest CP-even Higgs is SM Higgs lik^E Thus, we use the LEP II SM Higgs 
mass bound (rriH SM > 114.4 GeV) [26] as the lower bound on the lightest CP-even Higgs mass. The 
approximate one-loop expression of m^ is given by |27^l 

2 2 2~n 3G F mfr m 2 ~ A? f A? x i 
m^m|cos 2 2/3 + ^|-^ log -| + -± (l - ^) , (2) 



7T^ 



mi m~ 12m 



where m^ ~ m i R — m t L - ^ ne 1°S term in the bracket slowly increases as m^ increases. The second 
term proportional to (A t /mi) 2 is maximised when \A t /m^\ = The Higgs mass constraint 

requires large 7713, m^ 2 and also large \Aq\. 

As we aim to explore the small 7773 and m^ 2 region, we optimise Aq by taking the above two 
constraints into account at each (7713, m^ 2 ) point. Figure [1] (left) shows the allowed/excluded 
regions on the (7713 — Aq) plane at m^ 2 = 150 GeV. The pink and green regions are excluded by the 
tachyonic stop and the Higgs mass constraints, respectively. To obtain viable model points with 
small 7773, we take Aq = —m^ — 350 GeV throughout our analysis. 



4 We have checked gzzh/gzzH SM > 0.978 is hold in our parameter region, where the gzzH SM (dzzh) is the 
ZZHsm (ZZh) coupling. 

5 This equation is however not used in calculating low energy particle spectra. Instead, we use S0FTSUSY [28 v3.1.7 
program. 
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Figure [T] (right) shows the allowed/excluded region on the (7713 — m^ 2 ) plane. The pink and 
green regions are, again, excluded by the tachyonic stop and the Higgs mass constraints. The purple 
region is excluded by the chargino mass bound m~± > 94 GeV [29] from the direct search conducted 
by LEP II. The red (orange) region is excluded because the LSP is the lighter stop (stau). 

As we mentioned in the Introduction, FCNC and CPV observables do not place constraints 
on this parameter plane. As we adopt the large tuq (tuq = 1.5 TeV) and the moderate tan/3 
(tan/3 = 10), it is difficult to explain the anomaly of (g — 2) M [30] in our setup. We do not address 
this issue in this paper. 

Now we discuss the cross sections and branching ratios in our scenario. Figure [2] shows the 
total SUSY cross section and the cross sections for the leading production processes, gg and 
at yfs = 7 TeV pp collision. The values are obtained using the PR0SPIN0 [3H E2 EES EH E3 ESJ v2.1 
program. The next-to-leading order (NLO) corrections are taken into account. The total SUSY 
cross section is not large. It can be as large as O(10) pb only in the m^ 2 < 160 GeV region or the 
7713 < 700 GeV and m^ 2 < 260 GeV region. This is because the two dominant production processes, 
qg and qq, in the CMSSM (apart from the tuq ^> m 1 / 2 region) are highly suppressed due to the 
large mo value. The dominant process then turns out to be gg in the small m^ 2 region and t\t\ in 
the small 777,3 region. The other processes are negligible across the (7713 — m^) parameter plane. 

In this scenario, gluinos decay utterly into third generation squarks if the two-body decay mode 
g t\t or g — )► bib is open. Figure [3] shows Br(g —> t\t (£*£)) and Br(g — >> bib (bib)). The branching 
ratios are calculated using the SUSYHIT [37] vl.3 program. The g t\i mode dominates throughout 
the parameter space. Around 7713 c± 700 — 1200 GeV, those branching ratios abruptly become zero 
because the two-body decay modes become kinematically forbidden. In 7713 > 700 — 1200 GeV 
region, gluinos decay to charginos or neutralinos via three-body decays through off-shell squarks. 
The branching ratio of g — >> x/3/3, where gluinos decay into third generation quark pairs, together 
with weak gauginos, via three-body decay, is also shown in Figure [31 In 7713 < 1100 GeV region, the 
g x/3/3 mode has sizable branching ratios since off-shell stops and sbottoms are still lighter than 
first two generation squarks. The Br(g Xi 2^) an d Br(g Xi^b) rapidly decrease compared to 
Br(g —> Xi 2^) because of phase space suppression due to the top mass. 

3 ATLAS &-jet search and 0-lepton search 

To assess the impact of the direct SUSY searches on the MUSM scenario, we use two ATLAS 
search results: the ATLAS "6-jet" search with 35 pb -1 data and the ATLAS "0-lepton" search with 
165 pb -1 data. 

The cuts adopted in the ATLAS 6-jet search are shown in Table [TJ This search examines two 
signal regions: 0-lepton region and 1-lepton region. The 0-lepton region targets the g — >> bib mode. 
It requires at least one 6-tagged jet and adopts higher pt and E™ ss cuts. Events are discarded if 
they contain more than zero isolated leptons. On the other hand, the 1-lepton region targets the 
g —> tit mode. It requires at least one isolated lepton and 6-tagged jet assuming leptonic top decays. 
The pt and E^ 1SS cuts are mild compared to the 0-lepton signal region. This is because events that 
undergo the g — > tit modes contain a large number of final state particles. The pt of each final 
state particle can therefore not be large on average. In the MUSM scenario, as shown in Section 2, 
Br(g — >> tit (t*t)) + Br(g — >► bib (bib)) = 100% in the small 7713 region. It is therefore reasonable to 
expect this search has a good sensitivity to the MUSM scenario. 

The cuts used in the ATLAS 0-lepton search are shown in Table El This search defines three 
signal regions: 2-jets, 3-jets and 4-jets regions. All regions require very high pt jets and large E^ lss , 
since a larger number of SUSY events are available compared to the 6-jet search, due to 165 pb -1 



5 



gg cross section [pb] 




tit* cross section [pb] 



200 400 




600 800 1000 1200 1400 

m 3 [GeV] 



Total cross section [pb] 



200 400 600 800 1000 1200 1400 

ins [GeV] 




200 400 600 800 1000 1200 1400 

m 3 [GeV] 



Figure 2: Production cross section of supersymmetric particles. 
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Figure 3: Gluino branching ratio (only for the channels where the third generation quarks are involved). If 
777,3 is large and is small, the gluino cannot decay into two particles. In the last figure, all decay process 
containing t or b quarks are added. 
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Signal Region 


0-lepton 


1-lepton 


number of b-jets 


> 1 


> 1 


Ef iss [GeV] 


> 100 


> 80 


Leading jet px [GeV] 


> 120 


> 60 


Second jet p? [GeV] 


> 30 


> 30 


Third jet p T [GeV] 


> 30 




A0(jet,^ iss ) min 


> 0.4 


> 0.4 


Ef ss /m eS 


> 0.2 




m e g [GeV] 


> 600 


> 500 


m T [GeV] 




> 100 


SM Background 


19.6 ±6.9 


14.7 ±3.7 


Observed Events 


15 
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Table 1: The cut used to define two signal regions of ATLAS 6- jet analysis [23] . 



Signal Region 


> 2 jets 


> 3 jets 


> 4 jets 


E$ ISS [GeV] 


> 130 


> 130 


> 130 


Leading jet pr [GeV] 


> 130 


> 130 


> 130 


Second jet pr [GeV] 


> 40 


> 40 


> 40 


Third jet p T [GeV] 




> 40 


> 40 


Fourth jet pr [GeV] 






> 40 


A</>(jet,£™ ss ) min 


> 0.4 


> 0.4 


> 0.4 


Ef ss /m eS 


> 0.3 


> 0.25 


> 0.25 


m eS [GeV] 


> 1000 


> 1000 


> 1000 


SM Background 


12.1 ±2.8 


10.1 ±2.3 


7.3 ± 1.7 


Observed Events 


10 


8 
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Table 2: The cut used to define three signal regions of ATLAS 0-lepton analysis 0. 
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Figure 4: Reproduction of the ATLAS observed exclusion limit for validation of our method. Solid (dashed) 
lines are our (ATLAS') exclusion contours. 



data. Events are discarded if they contain more than zero isolated leptons. A very high effective 
mass cut (m e ff = J2jLi \Pt \ + ^t 1SS > 1000 GeV (for TV-jets region)) is also adopted. This search 
currently places the strongest limit on the CMSSM parameter space. 

(i) 

ATLAS provides the number of observed events n o ^ s that made it past cuts and the expected 

(i) (i) 

SM backgrounds n b together with their systematic error a h for each search signal region. Here i 

(i) 

represents the search signal regions. The are calculated by adding the uncorrelated background 
systematic and the jet energy scale systematic in quadrature. Those numbers are listed in Table [1] 
and [2] for each search signal region. 

(i) 

At each SUSY model point, the predicted number of signal events n s can be calculated using 

(i) (i) 

a Monte Carlo simulation. If one observes a statistically significant excess of the ny + n y b J from 
n^g 5 one can reject the SUSY model point at some confidence level. ATLAS presents the 95% CL 
exclusion regions in the CMSSM (mo — ^1/2) plane at the tan/3 = 40, Aq = 0, ji > slice for the 
6-jet search and at the tan/3 = 10, Ao = 0, ji > slice for the 0-lepton search. The dashed curves 
in Figure [H show the ATLAS' 95% CL exclusion contours. As can be seen, the mo < 400 GeV and 
mi/2 < 400 GeV region is excluded by the 0-lepton search. In the small mi/2 ^ 160 GeV region, mo 
is excluded up to 1000 GeV by both the 6-jet and 0-lepton searches. 

4 Monte Carlo simulation and its validation 

Before discussing the constraints using the ATLAS searches on the MUSM scenario, we summarise 
our setup for event and detector simulations. We generate 10 000 SUSY events at various model 
points using the HERWIG++ [38] v2.5.0 Monte Carlo program. The SUSY sample is scaled so that 
the corresponding luminosities are 35 pb -1 (165 pb -1 ) for the 6-jet search (0-lepton search). We use 
NLO SUSY cross sections obtained by using the PR0SPIN0 program. To simulate detector effects, 
we use the DELPHES [39J v. 1.9 program. The parameters defined in the DELPHES are tuned for the 
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ATLAS' analyses. For example, we use a R = 0.4 anti-fc^ algorithm for a jet reconstruction and 
assume the 6-tagging efficiencies to be 50%, 10% and 1% for 6-jet, c-jet and light flavour or gluon 

(i) 

jets. We analyse the SUSY sample and estimate the number of expected events n K s that survived 
the cuts which were used in the ATLAS searches. 

The ATLAS papers provide their estimate for the signal cross sections. Variations of the renor- 
malisation and factorisation scales by a factor of two results in uncertainties of 16%, 10%, 15%, 
27% and 30% for gg, q(f*\ qg, and b\b\ processes, respectively. The PDF uncertainties are 
expected to be 11 — 25%, 5%, 15% and 7 — 16% for gg, q(f*\ qg and i\t\ (bib*) processes, respec- 
tively. Those values depend on production processes and mass spectra for superparticles. We use 
the same numbers for the scale variations and the PDF uncertainties as the ATLAS b-jet searches 
but simply take the constant values. We take 25% and 11% for the PDF uncertainties on the gg 
and t\t\ (bib*) processes, respectively, across the SUSY parameter space. We then define the cross 
section error cr cross by adding those errors in quadrature. 

ATLAS does not supply the other systematic errors on the signal, such as jet energy scale 
uncertainty, 6-tagging uncertainty and luminosity uncertainty. To model these uncertainties, we use 

(i) . . (i) 

a single constant error a , for each search signal region. Our approximate a§ is then constructed as 

V a cross + ( a i^) 2 - We vary o^} and choose reasonable values so that our exclusion contours match 
the ATLAS' contours well. We choose 50% (50%) for the 6-jet (l)-lepton search signal region and 
10%, 30%, 30% for the 0-lepton 2-, 3-, 4-jets search signal regions, respectively. We find that the 

(i) (i) 

exclusion contours are not so sensitive to a , because a cross provides a sizeable contribution to ay . 

Given information, n^ s , n£\ and we can compute the exclusion p- value |. We 

follow Ref [10]. The expectation value for observed events is given as 

\® = n «Xl + 5 a *®)+r$\l + 6 b *®) t (3) 

where the impact of systematic variations is accounted for by the nuisance parameters 5 S and 5^. 
Using Poisson probability, the probability of observing n events is given by 

P(») = 4rf M,t e -><"(AC>) V||g+€ , 

WW J„«(-5,-l/4") ./max(-5,-l/„») "\ 

with the normalisation 

JV»= T d8 s f . d5 h e-^ 5 l\ (5) 

Jmax(-5,-l/cr^ ) ) Jmax(-5,-l/cr^ ) ) 

where we have truncated the Gaussian modelling of the systematic errors at 5a in computational 
practice. The lower edge of the integration is restricted to keep the signal and background con- 

(i) 

tributions independently non-negative. Finally, given observed events n o £ , the exclusion p- value, 
defined as the cumulative marginalised likelihood, is obtained as 

Pexcl(ni; ) s ) = ^P(n). (6) 

71=0 

The 95% CL exclusion region corresponds to p exc i < 0.05. 

6 There are various decent statistical methods, for example CLs, and those provide slightly different exclusion 
limits. Because the aim of this paper is not to investigate the dependence on the statistical methods, we simply use 
the frequent ist p- value method in our analysis. 
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The solid curves in Figure H] show the 95% CL exclusion contours obtained by the above proce- 
dure. The ATLAS' 95% CL exclusion contours are shown with dashed curves on the same figures. 
It can be seen that our exclusion contours reproduce the ATLAS' contours well. 

5 Constraint from ATLAS searches on MUSM scenario 

We are now ready to examine the constraints of the ATLAS searches on the MUSM scenario. 
We focus on the MUSM (7713 — m^ 2 ) parameter plane in the mo — 1.5 TeV, m# — 7713, Aq — 
7713 — 350 GeV, tan/3 = 10 slice, where the phenomenological allowed region has been found in 
Section 2. We divide the (7713, m^ 2 ) plane into grids with (75, 20) GeV intervals. At each grid 
model point, we calculate the p exc \ following the procedure described in Section 4 for each search 
signal region. 

Figure [5] shows the p exc i distributions for the 6-jet search in the 0-lepton (left) and 1-lepton 
(right) signal regions. The 95% CL exclusion contours are also shown in Figure [5l As can be seen, 
the 95% CL exclusion region is very limited. It is compacted only in the m^ 2 < 220 GeV region. 
We attribute this to the suppression of the gg cross section. In the m^ 2 > 220 GeV region, the gg 
cross section is less than lpb as shown in Figure El and the produced number of gg pairs may not 
be sufficient in the 35 pb -1 data. 

The excluded regions obtained from the 0-lepton and 1-lepton signal regions are complementary. 
The former signal region excludes the 7713 > 850 GeV and m^ 2 ^ 150 GeV region, whilst the latter 
excludes the 650 GeV < 7713 < 900 GeV and m^ 2 < 200 GeV region. The 0-lepton signal region 
adopts larger pt, E™ 1SS and m e ff cuts. It therefore prefers the gg events that undergo three-body 
gluino decays such as g 6&X^ 2 ) because the number of the final state particles is expected to 
be relatively small and the averaged pt of each particle can be sufficiently large. Indeed, gluinos 
decay to three-body decay modes only in the 7713 > 700 GeV region (See Figure El). In contrast, the 
cuts defined in the 1-lepton signal region are designed for the g t\i mode leading to long cascade 
decays and isolated leptons from leptonic top/stop decays. 

Figure [6] shows thep exc i distributions for the 0-lepton search in the 2-, 3- and 4-jets signal regions. 
As can be seen, the constraints from these search signal regions on the MUSM parameter space are 
rather weak despite the larger luminosity data with 165 pb -1 . There are no 95% CL exclusion 
regions found on the (7713 — m^ 2 ) plane in the 2- and 4-jets region, whilst the 4-jets region excludes 
the 1220 GeV < 7713 < 1350 GeV and m 1 / 2 ^ 130 GeV regior£|. This result is not surprising. The 

7 We have also applied the ATLAS 2010 (35 pb -1 ) 0-lepton analysis J4 to the MUSM scenario and checked that 
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Figure 6: p e xci-value distribution of 0-lepton search regions for MUSM scenario. The green curve corresponds 
to the 95% CL exclusion contour. 

0-lepton search requires rather large pt, E™ ss and ra e ff cuts. On the other hand, in the MUSM 
scenario, the pt of each particle can not be significantly large because of the small gluino mass and 
a large number of final state particles. 

To illuminate this observation, we show distributions of the £^ 1SS , leading jet pt and m e ff (2- 
jets) before cuts for representative model points for the CMSSM (blue histograms) and MUSM 
scenario (green histograms) in Figure [71 The model points we have chosen are mo = 350 GeV, 
m 1/2 = 250 GeV, A = 0, tan/3 = 10, ji > for the CMSSM and m 3 = 350 GeV, m 1/2 = 250 GeV, 
Aq = —700 GeV, tan/3 = 10, ji > 0, mo = 1.5 TeV for MUSM scenario. As can be seen, the 
distributions for the MUSM model point are much softer than those for the CMSSM point. The 
red arrows represent the cuts adopted in the 0-lepton search. We can see that the cuts remove the 
majority of events for the MUSM model point. Especially, the m e ff cut significantly reduces the 
number of signal events. To exclude the MUSM scenario, lower pt, E™ ss and m e ff cuts are crucial. 

6 Summary and conclusion 

Recent ATLAS and CMS SUSY searches have pushed up the exclusion limits significantly in the 
CMSSM (mo — ^1/2) plane. It may be getting more difficult to achieve naturalness within the 
framework of the CMSSM. We examined the impact of the recent direct SUSY searches on non- 
universal sfermion mass models focusing on naturalness. We considered non-universal sfermion mass 
models where the third generation sfermions involved in the 10-plet of SU(b) have a different soft 
mass 777,3 from the other sfermions' soft mass tuq. We focused on the parameter region where the 
parameters that are relevant to naturalness (777,3, ^1/2 •> m H u i aO are around the weak scale and the 
other dimensionful parameters have larger mass scales (mo = 1.5 TeV, \A$\ > 500 GeV). We applied 
the ATLAS 6-jet search and the latest 0-lepton search to various MUSM model points and identified 
the 95% CL exclusion region in the (777,3 — m i/2) plane. 

We found the constraints on the (777,3 — m i/2) plane from the ATLAS searches are rather weak. 
Those searches do not exclude the m^ 2 > 220 GeV region independently of 777,3. This result can 
be attributed to the following two features of the MUSM scenario. First, because of the large mo 
value, the cross sections of qg and qq processes are significantly reduced compared to the CMSSM 
(apart from the mo ^> m^ 2 region). Second, the main branching ratio of gluinos is typically 
g — > t\t —> xtbi — > x?^ + ^5 and events with such long cascade decays possess a large number of 
final state particles. Therefore, the pt of each particle can not, on average, be large. Most of such 
events fail to pass high pt cuts and large E™ ss and m e ff cuts. In this paper, we investigated the 
exclusion limit in the particular slice of the parameter space in the MUSM. The complete search 

there are no constraints on our (7713 — m 1 / 2 ) plane at the 95% CL. 
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Figure 7: Distributions of the £?™ 1SS , leading jet pt and m e ff for SUSY signals before event selection. 



across the whole parameter space is beyond the scope of this paper. However, as long as the 1st and 
2nd generation squarks are decoupled and the gluinos predominantly decay to the 3rd generation 
quarks, one can expect that the constraint is weaker than the CMSSM because of the reasons 
mentioned above. 

The non-universal sfermion mass models of this type can easily evade the current direct SUSY 
search constraints and keep naturalness. In order to exclude/discover these models, ordinary selec- 
tion cuts, based on high jets and large E™ ss and m e ff , are not efficient. The cuts based on the 
number of 6-jets, isolated leptons and the total number of final state particles may be preferable. 
We leave this study for future workJ§ 

8 During the publication process, ATLAS updated their analysis with the 1 fb _1 luminosity data collected in 2011. 
In the 0-lepton analysis [41] the exclusion region is substantially improved in the large mo region by raising the pr cut 
of the 4-jet signal region. However, one may naively expect that this update does not change our conclusion, because 
as we have seen in Figure[71 raising the pr cut significantly looses the efficiency for the MUSM scenario, although the 
detail study is required for the quantitative conclusion. For the 6-jet analysis, the constraint from the 1-lepton signal 
region on the simplified model, where gg and tit* productions are assumed and their decay chains are fixed, turned 
out to be weaker than for the 35 pb -1 result, since they observed an upward fluctuation in the data 02]. For the 
0-lepton signal region 03], the gluino mass limit for the simplified model, where gg and bib\ productions are assumed 
and their decay chains are fixed, is updated from 600 GeV to 700 GeV, by requiring more than or 1 6-tagged jets. 
This should extend the 95% exclusion region in the large wis region in Figure 5, although this region is less interesting 
compared to the small region where the fine-tuning problem is relaxed. The MUSM scenario may still provide an 
interesting possibility for the tension between naturalness and LHC constraint. Updating our analysis with the other 
LHC analyses is also our future work. 
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